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We study the interplay of the Kondo effect and spin-polarized tunneling in a class of systems 
exhibiting uniaxial magnetic anisotropy. Using the numerical renormalization group method we cal- 
culate the spectral functions and linear conductance in the Kondo regime. We show that the exchange 
coupling between conducting electrons and localized magnetic core generally leads to suppression of 
the Kondo effect. We also predict a nontrivial dependence of the tunnel magnetoresistance on the 
strength of exchange coupling and on the anisotropy constant. 
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Incorporating single atoms or molecules into nano- 
electronic devices is a very promising challenge, partic- 
ularly for information storage and processing technolo- 
gies [l|-|4( . Owing to recent achievements in experimental 
techniques, measuring transport through individual nat- 
ural and artificial (quantum dots) atoms and molecules 
has become feasible. Among various nanoscopic systems, 
magnetic atoms of spin S > 1/2 (like Fe, Co or Mn) 0- 
[f| and single- molecule magnets (SMMs), both exhibiting 
magnetic anisotropy 0, [J, have attracted much interest 
from both the fundamental as well as application points 
of view. In particular, it has been suggested that mag- 
netic states of SMMs can be controlled by spin-polarized 
current 0"H|, i.e. a sufficient current pulse can switch 
the magnetic moment between two low energy states 
(the model applies also to electric control of magnetic 
atoms). This has been recently proven experimentally 
in the case of magnetic adatoms [J| using spin-polarized 
STM technique. Furthermore, a SMM weakly coupled to 
two nonmagnetic metallic electrodes was shown to act as 
a spin filter [l^, while coupled to electrodes with differ- 
ent spin polarizations reveals features typical of a spin 
diode (l3|. In the strong coupling limit, on the other 
hand, the Kondo correlations play a significant role and 
anomalous features of transport become revealed at low 
temperatures. As shown experimentally, the magnetic 
anisotropy in such systems can be used to tune the Kondo 
effe ct Ifj , Il4| . Using break-junction technique, Parks et 



al. [14| were able to tune the anisotropy constant and 
therefore modify the energy spectrum responsible for the 
Kondo state, which in turn resulted in a crossover from 
the fully screened to underscreened Kondo effect. The 
Kondo phenomenon in transport through SMMs has been 
a subject of current interest, but the research so far was 
mainly focused on the interplay of the Kondo effect and 
quantum tunneling of SMM's spin [15| . 



To our knowledge, the physics of Kondo correlations 
in spin-polarized transport through SMMs and/or mag- 
netic adatoms is still rather unexplored. This problem is 
therefore addressed in the present Letter. We consider a 



situation when transport occurs via a local orbital of the 
system (orbital of a SMM, adatom, or a quantum dot), 
which is coupled to electrodes and additionally exchange- 
coupled to the corresponding magnetic core. The role of 
magnetic anisotropy in conductance and tunnel magne- 
toresistance (TMR) is also analyzed. We show, that the 
exchange coupling to the corresponding magnetic core 
generally suppresses the Kondo effect. Similarly, the mag- 
netic anisotropy has a significant impact on the conduc- 
tance and TMR in the Kondo regime. 

Model - We consider a generic model that includes es- 
sential features of such quantum objects like SMMs, mag- 
netic adatoms, and quantum dots exchange coupled to lo- 
cal magnetic moments, see Fig. [Ha). These systems are 
referred to as magnetic quantum dots (MQDs). We as- 
sume that MQD is coupled to ferromagnetic leads whose 
magnetizations can form either parallel (P) or antiparal- 
lel (AP) configuration, while MQD's magnetic easy axis 
is collinear with magnetic moments of the leads. We also 
assume that only a single orbital level (OL) of a MQD 
(lowest unoccupied molecular orbital of a SMM, atomic 
or quantum dot level) is directly coupled to the leads, 
while the corresponding magnetic core is coupled to the 
leads indirectly via exchange coupling to electrons in the 
local OL. The corresponding Hamiltonian reads 0, [lfl | 

Hmqd = -DS 2 Z + V e n a + U n tni - Js ■ S, (1) 

with S z denoting the zth component of the MQD's core 
spin operator S, and D being the uniaxial anisotropy con- 
stant of the MQD. The operator n a = c\c a . where c\(c a ) 
creates (annihilates) an electron of energy e in the OL, 
whereas U denotes the Coulomb energy of two electrons 
occupying the OL. Energy (and therefore occupation) of 
this level can be controlled by an external gate voltage. 
Finally, the last term accounts for exchange coupling be- 
tween the MQD's magnetic core and electrons in the OL, 
with s = (1/2) 4°V<r' c <r< and CT = (o- x ,<7 y ,o- z ) de- 
noting the Pauli matrices. 

Method - To accurately address the problem of trans- 
port through MQDs in the strong coupling regime, we 
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Figure 1. (Color online) (a) Schematic of a molecular quan- 
tum dot (magnetic molecule, magnetic adatom, or quantum 
dot exchange-coupled to a magnetic moment) coupled to two 
ferromagnetic electrodes, (b)-(e) Normalized spin-resolved or- 
bital level (OL) spectral functions a a (to) = TrT a A a (u) in the 
antiparallel (left) and parallel (right part of each figure) mag- 
netic configuration of electrodes for two different values of the 
exchange parameter J. The bottom curves (black) in each plot 
correspond to J = 0, with the solid line referring to a =f and 
the dashed one to a =\.. To increase readability, for J ^ 
the curves for a =\ [a =\) are shifted up and right by 0.2 
(0.4). The parameters are: S = 2, U = 1 meV, s/U = -0.5, 
T/U ~ 0.075, D/U « 1.7 • 10" 4 , and p = 0.5. 



use the Wilson's numerical renormalization group (NRG) 
method (lq . The NRG Hamiltonian of the full system 
can be then written as llTjl 



U = H M QD + Yl a V T */(np) [4/o<r + fL c a] 

+ _ [/na/n+lcr + fn+lafna] i (2) 

* — ^a.n— 

where /„ CT represents the nth site of the Wilson's 

chain (last term of H) and i„ denotes the hopping ma- 
trix element between neighboring sites of the chain. The 
second term of the NRG Hamiltonian stands for the cou- 
pling between the MQD and conduction electrons. Here, 
we use the flat band approximation, with the density of 
states p = 1/(22)), where 53 is the band half-width. The 
effect of ferromagnetic leads is then determined by the 
hybridization function r<j = T\^ a + Tro- , where r rcr is 
the coupling to the rth lead. When assuming left-right 
symmetry, the resultant coupling in the AP magnetic 
configuration does not depend on spin, P^m = T, and 
the system behaves then as being coupled to nonmagnetic 
leads. This is not the case in the P configuration where, 
rjj^jx = r(l±p), with p being the spin polarization of the 
leads and V = (T^ + TjJ/2. By solving Hamiltonian @ 



iteratively, we are able to determine static and dynamic 
quantities, basically at arbitrary energy fl7l | . Transport 
properties of a MQD are then determined from the spec- 
tral function of the OL, A a (uj) = — ilm{G^ (uj)}, where 
G^ (w) is the Fourier transform of the retarded Green's 
function, G*(t) = -i@{t){{c a {t), c\ (0)}). 

Spectral functions - The zero-temperature spin- 
resolved spectral functions A a (u>) of a singly occupied 
OL are shown in Figs. QJb)-(e) for the P and AP config- 
urations and for e = — U/2. To begin with, we note that 
in the AP configuration the effective coupling is the same 
for both spin orientations. This is, however, opposite to 
the P configuration where an effective spin splitting of the 
OL due to spin-dependent coupling to the leads (effective 
exchange field) occurs This exchange field depends 
on the position of the OL as fc eX ch ~ \n\e/(e + U)\. Con- 
sequently, the splitting will generally suppress the Kondo 
effect, except for the case of particle- hole symmetry point 
e = — U/2 (shown in Fig.[l| where fe C xch — ^ 0. Therefore, 
for to < Tk, where Tk is the Kondo temperature [ijj], the 
spectral functions in both configurations display a clear 
Kondo- Abrikosov-Suhl resonance at the Fermi level. 

When considering the effect of finite exchange interac- 
tion J, one should note that there are two competing in- 
teractions corresponding to the energy scales set by J and 
Tk- When \J\ < Tk, the system then tends to lower its 
energy by Tk due to formation of the many-body Kondo 
state. However, once \J\ > Tk, it becomes energetically 
more favorable for an electron in the OL to hybridize with 
the core spin 5* instead of free electrons in the leads. Thus, 
the resonant peak at the Fermi level will not develop in 
such a case. In Figs. QJb)-(e) we show the spectral func- 
tions corresponding to both ferromagnetic ( J > 0) and 
antiferromagnetic (J < 0) exchange couplings. It can be 
clearly seen that the Kondo effect becomes suppressed by 
the exchange coupling J, and this suppression is stronger 
for antiferromagnetic coupling. Actually, full suppression 
appears when J| > Tk- 

Linear conductance - The spectral functions deter- 
mine the linear conductance g a : 9^n\ = { e2 /h)(l — 
p 2 )nTA^^ (0) for the antiparallel configuration, and 

P /AP 

where A a (0) is the spectral function in the P/AP 
configuration for u> = 0. The linear conductance as a 
function of the OL energy is shown in Figs. Hta)-(f) 
for different values of J. First, for |J| <C Tk we re- 
cover results known for single- level quantum dots (20| . 
For the AP magnetic configuration one observes then an 
enhanced conductance due to the Kondo state when the 
orbital level is singly occupied. The conductance is then 
given by, fl^m = { e2 /h)(l — p 2 ), i.e. it is reduced by a 
polarization-dependent factor (1 — p 2 ) = 3/4, as com- 
pared to the conductance quantum, see Figs. [2ja,b) and 
also Figs. [3£a,b) for |J|/Tk <C 1. In the P configuration, 
in turn, the Kondo effect is suppressed due to effective 
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Figure 2. (Color online) (a)-(f) Dependence of the linear conductance g a on the OL energy e for indicated values of the exchange 
parameter J in the antiparallel (a)-(b) and parallel (c)-(f) magnetic configuration. Figures (g) and (h) show the corresponding 
TMR. Top panel represents the case of ferromagnetic (J > 0) coupling between electrons in the OL and MQD's magnetic core, 
whereas the bottom one corresponds to the antiferrornagnetic (J < 0) coupling. Parameters as in Fig. [T] 



spin splitting of the OL [18|, except for e = —U/2 where 
a sharp peak occurs and the Kondo effect is restored. 

When J / 0, the Kondo anomaly in conductance be- 
comes gradually suppressed with increasing |J|, and this 
suppression is faster for J < than for J > 0. This 
behavior is a consequence of a difference in quantum 
states taking part in formation of the Kondo state for 
J < and J > 0. First, the ground state energy Eqs 
of a singly-occupied bare MQD is lower for J < than 
for J > 0, £q F < Tgg, while the energies of virtual 
states corresponding to empty and doubly occupied OL 
are independent of J. Second, the singly occupied MQD's 
ground state for J < involves a superposition of both 
electronic spin states 'up' and 'down', which is not the 
case for J > 0. Consequently, the cotunneling processes 
driving the Kondo effect for J < are suppressed more 
effectively than for J > 0, and this behavior is clearly 
visible in the linear conductance, see Figs. Hta,b). When 
the Kondo peak becomes suppressed, only two main res- 
onances appear in the conductance, whose position de- 
pends on J: for J > they appear at e = JS/2 and 
e = —JS/2 — U, while for J < th eir p osition depends 
linearly on J and also weakly on D [10(. It is also inter- 
esting to note that in the P configuration the resonance 
peaks are almost exclusively due to spin-up conductance 
g P , see Figs. [Hc-f). This is related to the fact that spin- 
up electrons are the majority ones in both the left and 
right lead and thus tunneling of spin-up electrons is fa- 
vored, irrespective of J. The explicit dependence of the 
conductance on J is shown in Figs. [3Ja,b). 

Tunnel magnetoresistance - The effects associated 
with exchange coupling J are also pronounced in TMR, 



TMR = (g p - 5 AP )/.9 AP with 5 p / AP = £ CT .g^ /AP , see 
Figs. [2jg)-(h). For \J\ <C Tk, the conductance in the 
Kondo regime is generally larger in the AP configuration 
than that in the P one, g AP > g p , leading to negative 
TMR in major part of the Coulomb blockade regime, ex- 
cept for e « —U/2, when g AP < g p , and TMR is positive. 
When | J | increases, the Kondo peak in AP configuration 
becomes suppressed too, and positive TMR in the block- 
ade regime is restored. Since the suppression of g AP is 
more pronounced for J < 0, the corresponding TMR in 
the Kondo regime is larger for J < than for J > 0. 
Moreover, TMR for J < significantly exceeds the cor- 
responding Julliere's value, 2p 2 /(l — p 2 ) (= 2/3 for as- 
sumed parameters) [2l|. On the other hand, for empty 
or doubly occupied OL, where transport is mainly deter- 
mined by elastic cotunneling processes, one always finds 
g p > g AP with TMR approaching the Julliere's value. 

The explicit variation of the spin-dependent conduc- 
tance and the corresponding TMR with the parameter 
J is shown in Fig. [3] for different values of e. When 
e = -U/2 and | J\ -)• 0, TMR « p 2 /(l - p 2 ). In turn, as 
|J| grows to energies corresponding to the Kondo tem- 
perature, \J\ w Tk, TMR reaches a minimum, where it 
takes negative values. Further growth of \J\ above Tk is 
then accompanied by a significant increase in TMR, es- 
pecially for J < 0. In the Coulomb blockade regime with 
e 7^ —U/2, TMR is negative and constant for \ J\ < Tk, 
and starts increasing when \J\ > Tk to reach large posi- 
tive values for \ J\ ;§> Tk in the case of J < 0, see Fig.(3^b). 
In turn, at resonance, e = 0, TMR is positive and con- 
stant and only slightly decreases when | J| ~ T, where a 
local minimum develops due to the dependence of reso- 
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Figure 3. (Color online) Conductance g = ^ CT g a (a,b) and 
TMR (c,d) in the parallel (P) and anti-parallel (AP) configura- 
tions as a function of the exchange parameter J for indicated 
values of the OL energy e. The dashed vertical line indicates 
the Kondo temperature. Other parameters as in Fig. [T] 
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Figure 4. (Color online) Influence of the molecule's magnetic 
anisotropy D on TMR for the ferromagnetic (a) and antifer- 
romagnetic (b) exchange coupling J. Parameters are the same 
as in Fig. [T] with \J\ = 2T K . 

nance energies on J. 

The uniaxial anisotropy D modifies the energy spec- 
trum and electron states of MQD, and therefore has a 
significant influence on the Kondo state [ljj], which can 
be observed in the behavior of g and TMR. In Fig. [4] 
we show the e-dependence of TMR for different values 
of the anisotropy constant D. We note that the effects 
due to variation of D are mainly visible in the Kondo 
regime, while outside the Coulomb blockade TMR is 
rather independent of D. Moreover, this effect is more 
pronounced for J < than for J > 0. This is because 
magnetic anisotropy changes quantum states responsible 
for the Kondo effect, suppressing the difference in TMR 
for J < and J > 0, see Fig. g^b). 

Conclusions - Using numerical renormalization group 
method we have studied spin-dependent transport 
through a localized orbital level coupled directly to ferro- 
magnetic leads and exchange-coupled to a magnetic core. 
We have shown that there is a competition between in- 



teractions corresponding to two energy scales: the Kondo 
temperature Tk and exchange coupling J. The linear con- 
ductance becomes suppressed when |J| > Tk and this 
suppression is more pronounced in the case of antiferro- 
magnetic coupling J. Moreover, J also has a significant 
influence on TMR, which displays a nonmonotonic de- 
pendence on J with a minimum for \J\ « Tk and may be 
greatly enhanced when J < as compared to the case of 
J > 0. 
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